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Abstract 

Based on the formula for the number density of vacancies in a solid under the stress or tension, 
the model of grain boundary diffusion in crystalline solids is developed. We obtain the activation 
energy of grain boundary diffusion (dependent on the surface tension or the energy of the grain 
boundary) and also the distributions of vacancies and the diffusing species in the vicinity of the 
grain boundary. 



PACS numbers: 61.72.Bb, 66.30.Dn, 68.35.-p 



Recently, it was shown that sufficiently high pressures as well as mechanical stresses 
applied to a crystalline solid lead to the decrease in the energy of the vacancy formation 
and create, therefore, an additional amount of vacancies in the solid [1]. The last effect 
enhances self-diffusion in the crystal which is normally vacancy-mediated, at least in simple 
metals. Since large mechanical stresses are normally present in grain boundaries, these new 
results can elucidate the mechanisms of grain boundary diffusion which have remained so 
far unclear [2]. 

According to the thermodynamic equation [3] 

dE = TdS-pdV } (1) 

where E is the energy, T is the temperature, S is the entropy, p is the pressure, and V is 
the volume of a solid, the energy of a solid increases with pressure, so the pressure acts as 
the energy factor similarly to the temperature. Therefore, the number of vacancies in a solid 
increases both with temperature and with pressure. 

The thermodynamic consideration based on the Clausius- Clapeyron equation gives the 
number density n of vacancies in a solid in the form [1] 

n = (P /T) exp (-E v /T) = (n T /T) exp (-E v /T) , (2) 

where E v is the energy of the vacancy formation, P = n T is a constant, T can be put 
equal to the melting temperature of the solid at ambient pressure, and the constant no has 
an order of magnitude of the number density of atoms in the solid. Here the Boltzmann 
constant ks is included in the definition of the temperature T. 

The formula ^2} describes the thermal expansion of the solid. It should be taken into 
account that the dissolution of the vacancy gas in a solid causes the deformation of the 
crystalline lattice and changes the lattice parameters. 

The energy of the vacancy formation E v depends linearly on the pressure P (in the region 
of high pressures) as given by the formula 

E V = E - aP/n , (3) 

where a is a dimensionless constant, a ~ 18 for sufficiently high pressures. On the atomic 
scale, the pressure dependence of the energy of the vacancy formation in the equation is 
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produced by the strong atomic relaxation in a crystalline solid under high pressure. 

With increasing pressure, the number density of vacancies in a solid increases, according 
to the relation 

n = (n T /T) exp (- (E - aP/n ) /T) , (4) 

and, finally, the vacancies can condense, forming their own sub-lattice. Such is the explana- 
tion of the appearance of composite incommensurate structures in metals and some other 
elemental solids under high pressure [4-7]. 

Further increase of the number density of vacancies in a solid with increasing pressure 
leads to the melting of the solid under sufficiently high pressure (and fixed temperature). 
Such effect has been observed in sodium [6]. In general, such behavior is universal for solids, 
though the corresponding melting pressure is typically much larger than those for sodium. 

We assume that the melting of the crystalline solid occurs when the critical number 
density n c of vacancies is achieved. In view of the equation (j2j), it means that the ratio 
of the energy of the vacancy formation E v to the melting temperature T m of the solid is 
approximately constant, 

E v /T m « a. (5) 

The value of the constant a in the last relation can be determined from the empirical 
relation between the activation energy of self diffusion (which is approximately equal to the 
energy of vacancy formation) and the melting temperature of a solid [8]: 

E « 18T m , (6) 

so that a ~ 18. 

Substituting the expression (jHJ) in the relation (JHJ), we obtain 

(E - aP/n ) /T m « a. (7) 

The last equation gives the melting curve of the crystalline solid in the region of high 
pressures in the form 

T + P/no- E /a^T , (8) 
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where To is the melting temperature of the solid at ambient pressure. 

The constant no can be determined from the relation between the tensile strength a s and 
the melting temperature T m of a solid [1] 

n = cr s /T m . (9) 
The numerical value of this constant is no ~ 1.1 x 10 22 cm~ 3 [1]. 

Replacing in the relation (j3J) the pressure P by the absolute value of the stress or tension 
a = F/S, applied to a solid, where F is the applied force and S is the cross-section area of 
the solid in the plane perpendicular to the direction of the applied force, we can estimate 
the mean number density of vacancies in the solid under the stress or tension: 

(n) S (n T /T) exp (- (E - aa/n ) /T) . (10) 

The dissolution of the vacancy gas in a solid under the stress or tension is responsible 
for the low values of the elastic limit and the tensile strength of solids as compared with 
theoretical estimations not taking into account this process [9]. 

As indicated above, large mechanical stresses are normally present in grain boundaries. 
The absolute value o\, of the mechanical stress in the close vicinity of a grain boundary is 
given by the formula 

<Jb = lb/ro, (11) 

where 7& is the energy of the grain boundary and r is the radius of the atomic relaxation 
region (around a vacancy) which will be estimated below. 

According to the relation ffTOl) . the energy of the vacancy formation in the close vicinity 
of the grain boundary is given by the formula 

E b = E Q - aj b / (n r ) . (12) 

For the small values of misorientation angle 9 ^ 10 — 15 degrees, the energy of the 
dislocation structure contributes to the energy of the grain boundary [10]. However, for 
larger misorientation angles, the energy of the grain boundary is approximately constant 
and is determined by the surface tension 7 of the solid, 7^ = 7. 
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Due to the Einstein relation between the mobility of an atom, \i = v/F, where v is the 
velocity of the atom and F is the force acting on the atom, and the diffusion coefficient D 
[8]: 



H = v/F = D/T, 



(13) 



the speed of grain boundary motion v is proportional to the diffusion coefficient D± for 
self- diffusion in the direction perpendicular to the plane of a grain boundary. Therefore, 
the activation energy E of grain boundary motion is equal to the activation energy E± of 
self- diffusion across the grain boundary. The last activation energy is equal to the activation 
energy Ej, of grain boundary self-diffusion in the case of high-angle grain boundaries, and 
is approximately equal to the activation energy Eq of bulk self- diffusion for low-angle grain 
boundaries. Thus, there is a step of the activation energy for grain boundary motion at 
some critical value 9 C of the misorientation angle (9 C = 10 — 15°, as indicated above). Such 
a step of the activation energy for grain boundary motion has been observed experimentally 
in high-purity aluminium, the critical value of the misorientation angle being in this case 

e c = 13.6° [li]. 

The driving force for grain boundary motion is provided by the distribution of mechanical 
stresses in a crystalline solid [12]. 

Assuming that the free surface of a crystalline solid is formed by the plane of vacancies, 
we can estimate the surface tension of the solid as follows 



where a = n = 0.45nm has an order of magnitude of the lattice spacing a, and (3 is a 
dimensionless constant which has an order of unity. For hard metals such as Al, Zr, Nb, Fe, 
Pt, f3 = 0.8. In the case of mild metals, (3 is normally smaller, e.g. for Rb and Sr, /? = 1/4. 

Substituting the estimation ( fl4l) for the energy of the grain boundary 75 = 7 in the 
equation (1121) . we find 



Due to the atomic relaxation and thermal motion of atoms, the migration barriers are 
small [2,13], and the activation energy of self-diffusion is approximately equal to the energy 



7 = (3n E a , 



(14) 




(15) 
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of the vacancy formation. The analysis of experimental data on the activation energy of 
grain boundary self-diffusion gives an empirical relation [14] 

E b « 9T m « E /2. (16) 

From equations (Tl5!) and (TT6l) . we find the estimation of the radius of the atomic relaxation 
region, 

r ~ 2(3aa = aa , (17) 
since j3 has an order of unity. The radius of the atomic relaxation region has an order of 

— 1/3 

r = 18n m 8nm. This value is comparable with the diameters of tracks produced by 
high energy ions in metals [15-17]. The grain boundary diffusion width 5 [14] is smaller than 
the radius of the atomic relaxation region due to the non-uniform distribution of vacancies 
inside the atomic relaxation region in the grain boundary. 

If we assume that the mechanical stress a decreases linearly with the distance x from the 
plane of the grain boundary, 

o = a (1 - kx) , (18) 

where o"o is the stress at the boundary of the atomic relaxation region with the width r in 
the grain boundary (this value is smaller than o b = 7/r = (1/2) noT m and has an order of 
magnitude o"o = (1/2) n T), then the equation (flOj) gives the distribution of vacancies in the 
vicinity of the grain boundary in the form 

n = {jiqTq/T) exp (— (E — aa (1 — kx) /n ) /T) = n b exp (—aa kx/ (n T)) , (19) 

where rif, is the number density of vacancies at the boundary of the atomic relaxation region. 

Due to the trapping by vacancies [18], the distribution of the concentration c of the 
diffusing species in the vicinity of the grain boundary follows the same law: 

c b exp{-x/l), (20) 

where q, is the concentration of the diffusing species at the boundary of the relaxation region, 
and the scale / is given by the formula 
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/ = n T/ (aa k) . (21) 

Here k has an order of magnitude of 1/d, d being the size of the grain, so that I = d/a. The 
penetration profiles described by the equation (1201) have been indeed observed experimentally 
in the case of grain boundary diffusion in metals [8, 18], the measured penetration depth / 
having an order of a few micrometers [8]. 

To summerize, we obtained the dependence of the activation energy of grain boundary 
self-diffusion on the energy of the grain boundary, the estimation of the surface tension of a 
solid and of the energy of the grain boundary, and the width of the atomic relaxation region 
in the grain boundary (or the radius of the atomic relaxation region around a vacancy). We 
obtained further the distributions of vacancies and the diffusing species in the vicinity of 
the grain boundary. The obtained radius of the atomic relaxation region is consistent with 
the diameters of tracks produced by high energy ions in metals. 
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